Piling and avalanches of magnetized particles 
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We performed computer simulations based on a two-dimensional Distinct Element Method to 
study granular systems of magnetized spherical particles. We measured the angle of repose and the 
surface roughness of particle piles, and we studied the effect of magnetization on avalanching. We 
report linear dependence of both angle of repose and surface roughness on the ratio / of the magnetic 
dipole interaction and the gravitational force (interparticle force ratio). There is a difference in 
avalanche formation at small and at large interparticle force ratios. The transition is at f c m 7. For 
/ < f c the particles forming the avalanches leave the system in a quasi-continuous granular flow 
(granular regime), while for / > f c the avalanches are formed by long particle clusters (correlated 
regime). The transition is not sharp. We give plausible estimates for f c based on stability criteria. 
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I. INTRODUCTION 

The dipole interaction between magnetized particles 
can be viewed as an anisotropic adhesion force. Because 
its strength can be easily manipulated by the strength 
of the magnetizing field, magnetized particles have been 
recently proposed pj, |2j to get insight into the transition 
from noncohesive to cohesive grains. Previously adhesion 
effects have mainly been studied in form of moisture- 
induced changes, significant for industrial processes in 
fields as pharmaceuticals, agriculture, and constructions. 

Some time ago Hornbaker et al. Q addressed the 
question how sand castles stand. They stated that al- 
ready small quantities of wetting liquid can dramatically 
change the properties of granular media, leading to large 
increase in the angle of repose and correlation in grain 
motion. Theoretical studies on the angle of repose based 
on stability criteria have been done by Albert et al. £|. 
They theoretically determined the dependence of the an- 
gle of repose on cohesive forces, and applied the results 
to wet granular material. 

Experimental studies of Tegzes et al. on angle of repose 
using the draining crater method p| and on avalanches 
using a rotating drum apparatus 6] identify three distinct 
regimes as the liquid content is increased: a granular 
regime in which the grains move individually, a correlated 
regime in which the grains move in correlated clusters, 
and a plastic regime in which the grains flow coherently 

Experiments of Quintanilla et al. pfl using the rotat- 
ing drum apparatus address the question of self-organized 
critical behavior in avalanches of slightly cohesive pow- 
ders. Their results show that avalanche sizes do not fol- 
low a power-law distribution, however, they scale with 
powder cohesiveness. Samadani et al. §| studied the 



effect of interstitial fluid on the angle of repose and the 
segregation of granular matter poured into a quasi-two- 
dimensional silo. 

To study the transition from noncohesive to cohesive 
behavior, Forsyth et al. 0,0, adopting the widely sug- 
gested idea that competition between the interparticle 
forces and the inertial forces determines the behavior of 
cohesive granular materials, suggested a method based 
on magnetized particles. The particles placed in an ex- 
ternal magnetic field become magnetized, all having the 
same magnetic orientation parallel to the field. Vary- 
ing the strength of the field allows to continuously vary 
the resulting interparticle magnetic force. Using nonmag- 
netic perspex walls the particle-wall interaction remains 
unchanged relative to the noncohesive state. Using par- 
ticles under same packing conditions it is ensured that 
the initial conditions are as uniform as possible. 

We carried out computer simulations on a system cor- 
responding to the experiments of Forsyth et al. Q, Q 
and studied the angle of repose, the surface roughness, 
and the effect of magnetization on avalanching in two- 
dimensional particle piles. 

The magnetic interaction of magnetized grains is 
highly anisotropic, and the fixed external field introduces 
even more anisotropy as the grains are aligned to the 
field. A similar experimental setup 9], but with parti- 
cles carrying a remanent magnetization in the absence 
of an external magnetic field, would partly diminish the 
mentioned anisotropy, however in this case the magneti- 
zations and the interparticle forces are not as well defined 
as in the experiments of Forsyth et al. 0, • 

Because of the strong anisotropy and the longer inter- 
action range one can expect differences between the re- 
sults on magnetized particles and wet granular systems, 
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however the basic effects should be the same. 



II. SIMULATION METHOD 

We performed computer simulations based on a two- 
dimensional Distinct Element Method (DEM) ^(} (f° r a 
review see 0, 0, and references therein) to study 
granular systems of magnetized spherical particles. The 
particles are magnetized by a constant external field, all 
having the same magnetic orientation parallel to the field. 
The magnetization is modeled with dipoles. We neglect 
any coupling between the magnetic orientation and par- 
ticle rotation (i.e., the particles can rotate freely, while 
their magnetic dipole is fixed). 

For characterizing the strength of the interparticlc 
force, we introduce a dimensionless quantity defined by 
the ratio of the maximum magnetic interparticle force at 
contact and the gravitational force. 

The magnetic force acting on a dipole m 2 situated at 
distance r 2 i from a dipole mi, along the direction n 2 i is 
given by 



llq 3 r . . , . 

F21 = -. 4- (n 2 im 2 ) mi + (n 2 imi) m 2 

4tt r\ x 

-5 (n 2 imi) (n 2 im 2 ) n 2 i + (mim 2 ) n 2 i ] . (1) 

For identical hard spherical particles of diameter D 
and magnetic dipole S, the largest possible dipole-dipole 
magnetic force is 
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which corresponds to a head-to-tail configuration, the 
dipoles having the same orientation. 

Wc define the magnetic interparticle force ratio as 



/ = F m /F, 
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where F g = mg is the gravitational force (m denotes the 
particle mass and g is the gravitational acceleration). 

Considering mass density p and magnetization M, we 
have m = pV, S = MV, V = irD 3 /6, and thus 
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Assuming some / interparticle force ratio, from the 
previous equation the corresponding magnetization can 
be calculated as 
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In our simulations we used p — 7.5 g/cm 3 (which cor- 
responds approximately to the mass density of steel), 
g = 9.8 m/s 2 , and interparticle force ratio / < 25. 

The diameter of the spherical particles was taken from 
the 0.7 — 0.9 mm interval, with a Gauss-like distribu- 
tion having the mean of 0.8 mm. The vast majority of 
particles had diameters very close to D = 0.8 mm. The 
Gauss-like distribution is given by the average of 4 inde- 
pendent uniformly distributed random variables in the 
mentioned interval. This slight polydispersity, resem- 
bling real experimental setups, is used to avoid effects 
originated from symmetries of monodisperse systems. 

The long range magnetic interaction is taken in con- 
sideration within a reasonable cutoff distance as a dipole- 
dipole interaction. We choose the magnetic interaction 
cutoff to 6.25D (where D is the average particle diame- 
ter). As shown in a previous study 0], 5D value already 
gives a reasonable magnetic interaction cutoff in two- 
dimensional dipolar hard sphere systems regarding the 
local ordering. The angle of repose, the surface rough- 
ness, and the particle avalanches depend crucially on lo- 
cal orderings inside the pile, as noted for example by 
Altshuler et al. in 15] . The used cutoff keeps the char- 
acter of local orderings and changes the magnetic energy 
per particle by less than 5% [l4|. 

We calculate the collision interaction of particles using 
the Hertz contact model |l6j with appropriate damping 
[l7|. We implement Coulomb sliding friction for large 
relative translational velocities and for numerical stabil- 
ity viscous friction for small velocities, with continuous 
transition between the two, controlled by a (large) vis- 
cous friction coefficient. We do not use any static or 
rolling friction model. A grid based method is used to 
identify neighboring (and potentially colliding) particle 
pairs. 

The parameters of the Hertz contact model were cho- 
sen such that they correspond to Young modulus of ap- 
proximately 0.015 GPa and restitution coefficient of ap- 
proximately 0.86. These are characteristic values for hard 
rubber elastomers (used for example in constructing golf 
ball covers). The Young modulus is orders of magnitudes 
smaller than that of steel, a choice enabling realistic CPU 
times with the DEM method. The particle-particle and 
the particle- wall sliding friction coefficient was 0.5 and 0.7 
(characteristic for steel-steel and steel-perspex friction) . 

The translational motion of particles is integrated 
based on Newton's equation using Ver let's leap-frog 
method. The rotational state of particles is integrated 
with Euler's method. The integration time step was 5 ps. 
With the used elastic parameters a good lower estimate 
for collision times is 170 ps. In such conditions, the used 
integration time step gave good numerical stability and 
also fairly good response time on PCs with 1.8 GHz CPUs 
available at the time of writing. 

The simulation setup can be seen in Fig. ^ The ex- 
ternal magnetic field is vertical. The particles are added 
one by one with constant rate along vertical trajectories 
at small (maximum one particle diameter) random dis- 
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FIG. 1: Simulation setup. The particles are introduced with 
constant rate one by one at small random distances from the 
left wall. The particles can leave the system on the right side. 
The surface angle (i.e. the angle of repose) is measured by 
fitting a straight line over the positions of the surface parti- 
cles (marked with black). The figure also shows the normal 
contact forces. The thickness of the lines connecting the cen- 
ters of the particles in contact is proportional to the normal 
contact force. The sample corresponds to / = 6 interparticle 
force ratio. For color figure see the electronic version. 

tance from the left wall. They either reach the pile with 
a given velocity (i.e. they are fired into the pile), or their 
impact velocity is set to zero (i.e. they are placed gently 
on the pile). 

The system's bottom wall is sticky. Any particle touch- 
ing the bottom wall sticks to the wall. This builds up a 
random base (see the experimental setup used by Alt- 
shuler et al. described in ^j|). The particles can leave 
the system on the right side. The particles are removed 
from the simulation when their distance from the bottom- 
right corner is larger than the magnetic interaction cutoff 
distance. The system's bottom is 51.251? length. Only 
this system size was used, finite size effects were not stud- 
ied. 

The surface particles (marked with black in Fig. ^) 
are identified with the weighted alpha shape algorithm 
[18L Il9| . Alpha shapes are generalizations of the con- 
vex hull and can be used for shape reconstruction from 
a dense unorganized set of data points. The weighted 
alpha shapes are extensions of this kind of shape recon- 
struction to a set of spheres (as in our case). We used 
the implementation included in the Computational Ge- 
ometry Algorithms Library poj . The algorithm's alpha 
parameter was set to the square of the mean particle size. 
This gave satisfactory results. The surface angle (i.e. the 
angle of repose) is measured by fitting a straight line over 
the positions of the surface particles. The surface rough- 
ness is given by the standard deviation of the surface 
points from the fitted line. 

As part of our investigations, with a special side wall 
model, we also simulated the effect of front and back walls 
in a Hele-Shaw cell geometry encountered in experimen- 



tal studies. We took into consideration the frictional in- 
teraction with side walls by summing the magnitude of 
normal forces acting on one particle, directing a well de- 
fined percentage of this pressure on the walls, and deriv- 
ing a frictional force using the already mentioned friction 
model. The percentage of the total force directed on the 
side walls was a parameter of our simulations. 

We performed three sets of simulations: (a) the parti- 
cles were fired into the pile, (b) the particles were placed 
gently on the pile, and (c) the particles were fired into 
the pile, while 4% of the internal pressure was directed 
on the front and back walls (see side wall model) . In both 
(a) and (c) the particles reached the pile with 0.5 m/s im- 
pact velocity, which corresponds to approximately 16D 
dropping height. In all three simulation sets we executed 
runs at different interparticle force ratios. In each run 
we started with an empty system, and introducing 12000 
particles, one particle every 3000 integration steps, we 
numerically integrated the system for 3 minutes (simu- 
lated time). 

In the first part of the process the number of particles 
in the system increased monotonically. After a pile was 
built, avalanches started, which in a pulsating manner 
moved particles out of the system. In this way the num- 
ber of particles began to oscillate around some well de- 
fined value. In this latter part we identified the surface 
particles every 500 integration steps, and we measured 
the slope of the fitted surface line and the standard de- 
viation of surface points from this line. The average of 
this quantities over the simulated time gave the measured 
angle of repose and surface roughness. 

We also measured the avalanche sizes and avalanche 
durations. This can be done in many different ways. We 
define an avalanche by a number of individual events on 
the time scale of the integration steps (i.e. the small- 
est simulation time step) in which (at least) one particle 
leaves the system, and the time between two consecutive 
events is smaller than a well defined value. We take this 
value equal to the time corresponding to 3000 integration 
time steps. Our choice is based on the system's observed 
dynamic time scale, and the fact that one new particle 
is introduced (i.e. the system is perturbed) every 3000 
integration time steps, and thus on a larger time scale 
there are surely uncorrelated events. 



III. SIMULATION RESULTS 
A. Angle of repose and surface roughness 

In all cases both the angle of repose and the surface 
roughness (in the examined domain) exhibit a linear de- 
pendence on the interparticle force ratio (see Fig. I2J • 

The angle of repose in case (a) and (b) increases by 
approximately 0.5 degree per unit change of interparticle 
force ratio (see upper panel in Fig. This is in good 
accordance with the experimental results of Forsyth et al. 
0, , however the angle of repose at zero magnetization 
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FIG. 2: Angle of repose (upper panel) and surface rough- 
ness (lower panel) at different magnetic interparticle force ra- 
tios. The angle of repose is measured in degrees. The surface 
roughness is measured in (average) particle diameters. We 
executed three set of simulations. In (a) and (c) the parti- 
cles were fired into the pile, while in (b) the particles were 
placed gently on the pile. In (c) an artificial side wall effect 
was switched on (see text for details). 



FIG. 3: Distribution of particle avalanche sizes (upper panel) 
and avalanche durations (lower panel) at zero magnetization. 
We examined three different simulation setups (see text for 
details) . Firing the particles into the pile (i.e. dropping them 
from a given height) or placing them gently, and switching on 
or off the side wall effect gave no qualitative difference. Over 
our simulation data (for both avalanche sizes and avalanche 
durations) we could fit a stretched exponential with 7 = 0.43. 



in our case is about 10 degrees smaller. This can be 
the result of the missing side wall effect (see for example 
|2lj)- and the missing static and rolling friction (see for 
example jUii). 

At zero magnetization the average surface roughness is 
about 0.7 particle diameters, and in all cases increases by 
approximately 0.12 particle diameters per unit change of 
interparticle force ratio (see lower panel in Fig. 01 • 

As a consequence of our side wall model, the angle of 
repose in case (c) is about 8 degrees higher at zero mag- 
netization in agreement with experimentally observed ef- 
fects of front and back walls in Hele-Shaw cells. How- 
ever, the way we model the side walls leads to a stronger 
increase of the angle of repose with / than in the exper- 
iments of Forsyth et al. 0, Q (see inset of upper panel 
in Fig. |2| . The side wall effect does not influence the 
surface roughness (see inset of lower panel in Fig. EI> • 



B. Particle avalanches at zero magnetization 

We carefully examined the distribution of particle 
avalanches at zero magnetization in all three simulation 
sets. We also executed extra runs introducing a total 
of 144000 particles, one particle every 3000 integration 
steps, integrating for a total of 36 minutes (simulated 
time). It must be noted that simulating at zero magne- 
tization (i.e. without magnetization) is about 10 times 
faster, because only the (short-range) collision interac- 
tion must be calculated. This permitted longer simula- 
tion times. The corresponding particle avalanche size and 
duration statistics can be seen in Fig. 

Firing the particles into the pile (i.e. dropping them 
from a given height) or placing them gently, and switch- 
ing on or off the side wall effect gave no qualitative differ- 
ence. Over both the avalanche size and duration distri- 
bution data we could fit stretched exponentials of form 
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P(x)=P exp[-(x/x y] (6) 

with 7 = 0.43. 

Our 7 value is in good agreement with the experimental 
results of I24J and is the same value as the one found by 
Frette et al. |23 f° r piles of rice with small anisotropy. By 
contrast recent work of Costello et al. |2(J presents care- 
fully collected and detailed experimental results on piling 
of uniform spherical glass beads, which show a power law 
behavior with an exponential cutoff. Costello et al. argue 
that the exponential cutoff depends on the height from 
which the particles are dropped and probably also on co- 
hesion forces. It was not our intention in this study to 
contribute to the clarification of this point. 



C. Effect of magnetization on particle avalanching 

We analyzed the effect of magnetization on particle 
avalanching, and we found that there is a difference in 
avalanche formation at small and at large interparticle 
force ratios (see avalanche movies |23). We identified a 
granular and a correlated regime. The transition between 
the two regimes is not sharp. Similar regimes were iden- 
tified experimentally by Tegzes et al. [3, in case of 
wet granular materials. At high liquid content (i.e. large 
interparticle force ratio) they could also identify a third 
plastic regime. 

Studying the recordings from our simulations |27| . 
it can be observed that for small magnetizations the 
avalanches are formed by small vertical chains follow- 
ing each other at short times, giving the impression of 
a quasi-continuous flow (granular regime). As the par- 
ticle magnetization increases, at f c ~ 7 the previously 
quasi-continuous flow is replaced by individual narrow 
and long particle clusters falling at the system's bound- 
ary (correlated regime). 

There are two elementary processes characteristic for 
the dynamics of avalanche formation: peeling and split- 
ting (see Fig. BJ. Near the foot of the pile the particles 
are arranged approximately in a triangular lattice with a 
horizontal base. By contrast magnetic interactions would 
favor a triangular lattice with a vertical base. Clusters 
of particles near the free surface of the pile can peel off 
the pile by rotating 30 degrees into this favorable con- 
figuration (see part (i) of Fig. 0J). A favorably oriented 
domain of particles on a triangular lattice consists of ver- 
tical chains shifted by half a particle diameter from one 
chain to the next. This domain may be stable or may 
disintegrate into smaller domains or chains, depending 
on whether or not the magnetic interactions are strong 
enough to prevent the dilation, which is necessary to al- 
low for relative motion of chains within the domain (see 
part (ii) of Fig. 0). 

Combinations of these two processes determine the size 
of the outflowing clusters, which typically consist of v 



(i) 



(ii) 




FIG. 4: Elementary processes characteristic for the dynam- 
ics of avalanche formation: (i) peeling and (ii) splitting. A 
particle cluster composed of parallel chains can peel off the 
pile by rotating to a more favorable orientation, and can dis- 
integrate into clusters composed of smaller number of chains. 
Combinations of these two processes determine the size of the 
outflowing clusters. 

parallel vertical chains of length N. Both N and v de- 
pend on /. Taking into account only the magnetic inter- 
action between nearest neighbors one can give plausible 
estimates for this dependence. 

We compare the gain in magnetic energy and the loss 
in gravitational energy when a cluster of particle chains of 
length N rotates by an angle a into a position where the 
chains are aligned with the magnetic field. The magnetic 
moment of each particle is always aligned with the mag- 
netic field, so that only the dipolar interaction between 
the grains matters. We take into consideration only intra- 
chain interactions, neglecting the interchain interactions. 
The dipolc-dipole-interaction potential between two par- 
ticles whose center of mass connection is tilted by an 
angle a with respect to the magnetic field is 

^ = ^5 (l-3eos 2 «). (7) 

Hence rotating them into the field direction (a = 0) low- 
ers the energy per particle in the cluster by an amount 
proportional to 

A£ m = £^J (l-cos 2 a). (8) 

The center of mass is lifted by this rotation. Therefore 
the gravitational potential energy per particle in the clus- 
ter increases by 

AE g = N (1-cosa). (9) 
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Setting a = 30 degree and using the definition of / given 
in Eq. Q, Ai? m = Ai? ff shows that clusters up to chain 
length 

(1 + V3/2) / (10) 

can peel off. 

In order to check to what thickness v max such a cluster 
is stable with respect to splitting into subclusters with 
less chains, we compare the magnetic energy loss per unit 
chain length with the gravitational energy gained, when 
a subcluster of v chains moves down by half a particle 
diameter relative to the rest of the cluster. In this case 
we have 

Un S 2 , . „ mqD 

AEm = £ 2D*' ^ AE ° = 2 V - 

This shows that the dipole-dipole-interaction can only 
prevent cluster splitting, if the cluster consists of less than 
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FIG. 5: Dependence of average avalanche flow on interparti- 
cle force ratio in granular regime. We examined three different 
simulation setups (see text for details). Firing the particles 
into the pile (i.e. dropping them from a given height) or plac- 
ing them gently, and switching on or off the side wall effect 
gave no qualitative difference. The avalanche flow is measured 
in particles per second. 



v max ~ g / ( 12 ) 

chains. 

Based on the above results we can discuss the pro- 
cess of avalanche formation. Already at small magnetiza- 
tions the surface roughness allows for coherent rotation of 
larger clusters up to chain length N max . For / < f c w 6 
clusters consisting of more than one chain of particles 
can easily dilate and will disintegrate into isolated chains, 
forming a quasi-continuous flow. For f > f c clusters con- 
sisting of up to v max > 1 chains can fall. These results 
are close to the observations made on the simulations: 
The transition between the two regimes was observed at 

The difference in avalanche formation in different 
regimes can be clearly observed in avalanche duration 
to avalanche size relation. In all three simulation sets, 
at given interparticle force ratios, for each avalanche size 
we collected the measured avalanche durations and cal- 
culated the corresponding average avalanche durations. 
We also examined the avalanche size and duration dis- 
tributions. Our results are summarized in the next two 
subsections. 



D. The granular regime 

In the granular regime the avalanche sizes are propor- 
tional to the corresponding average avalanche durations. 
The proportionality factor (i.e. the ratio of avalanche 
sizes and average avalanche durations) defines an aver- 
age avalanche flow, which increases linearly with / (see 
Fig. This linear dependence is explained by the fact 
that the avalanches in granular regime are formed by 



quasi-continuous flows of small particle chains, and the 
height of these chains increases linearly with /. 

Independent of the setup, at zero magnetization the 
average avalanche flow is approximately 182 particles per 
second, and increases with approximately 74 particles per 
second per unit change of interparticle force ratio (see 
Fig.0. 

The avalanche size distribution at interparticle force 
ratios 1 < / < 7 can be scaled together reasonably well 
(see Fig. |SJ using the ansatz 



^(s,/) = r 1 Q(.s//), (13) 

where s denotes avalanche sizes, P(s, /) is the probability 
associated with an avalanche size, and Q(-) is a function 
with integral 1 on the [0, +00) interval. 

Based on the avalanche size distributions (see Fig.EJ, 
we argue that the magnetic cohesion introduces a well- 
defined characteristic size in particle avalanches. From 
the scaling property, we conclude that, the characteris- 
tic avalanche size increases linearly with the interparticle 
force ratio. Qualitatively similar results were found in 
the experiments by Szalmas et al. 0. 

As both the characteristic avalanche size and the aver- 
age avalanche flow increase linearly with /, the charac- 
teristic average avalanche duration (equal to the ratio of 
the former two) in leading order is independent of /. The 
dependence of the avalanche duration distribution on / 
is contained in higher order corrections which could not 
be captured by our simulations. 
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FIG. 6: Scaled avalanche size distribution in granular regime. 
A well-defined characteristic size can be observed. We exam- 
ined three different simulation setups (see text for details). 
The avalanche size distribution at interparticle force ratios 
1 < / < 7 are scaled together using the ansatz P(s,f) 
f Q(s/ /), where s denotes avalanche sizes. From the scal- 
ing property, we conclude that, the characteristic avalanche 
size increases linearly with /. 



FIG. 7: Scaled avalanche duration distribution in corre- 
lated regime. A well-defined characteristic duration can be 
observed. We examined two simulation setups (see text for 
details). The avalanche duration distribution at interparticle 
force ratios 7 < / < 25 are scaled together using the ansatz 
P[ T :f) = f~ 1 ^ 2 Q( T / f 1 ^ 2 )) where r denotes avalanche dura- 
tions. From the scaling property, we conclude that, the square 
of characteristic avalanche duration increases linearly with /. 



E. The correlated regime 

In correlated regime the avalanche durations are given 
by the free fall of long particle clusters, and accordingly 
the square of the measured avalanche durations are pro- 
portional to the length of the clusters, and thus propor- 
tional to /. At the same time, as it was already men- 
tioned, the width of the falling particle clusters is small, 
and thus the avalanche sizes in leading order are propor- 
tional to the cluster length. In consequence the avalanche 
sizes are proportional to the square of the corresponding 
average avalanche durations, contrary to the linear de- 
pendence found in granular regime. 

As it was already mentioned, the width of the falling 
particle clusters can take arbitrary values up to some well 
defined maximum. In other words, the particle clusters 
can have different number of layers while having the same 
length. This introduces large fluctuations in avalanche 
sizes. These fluctuations are proportional to the cluster 
length, and thus proportional to /. 

The avalanche duration distributions at interparticle 
force ratios 7 < / < 25 can be scaled together reasonably 
well (see Fig. UJ using the ansatz 

P(r,f) = r l/2 Q(T/f^), (14) 

where r denotes avalanche durations, P(r, /) is the prob- 
ability associated with an avalanche duration, and Q(-) 
is a function with integral 1 on the [0, +oo) interval. 

A well-defined characteristic duration can be observed 
in Fig. [7| Based on the scaling property, the square of 
the characteristic avalanche duration increases linearly 
with the interparticle force ratio. As the square of the 



avalanche duration is proportional to a corresponding 
cluster length, and a cluster length in leading order de- 
fines an avalanche size, we can conclude that there is also 
a mean characteristic avalanche size, which consequently 
is proportional to /. We could not find this avalanche 
size explicitly in the avalanche size distribution because 
of the large fluctuations of the same order of magnitude 
as the mean value. 



IV. CONCLUSIONS AND DISCUSSION 

We studied basic effects in two-dimensional granular 
piles formed by magnetized particles in simulations sim- 
ilar to the experiments of Forsyth et al. Q, ■ We mea- 
sured the angle of repose, the surface roughness, and the 
effect of magnetization on particle avalanching. As a gen- 
eral result we can mention that dropping the particles 
from a given height or placing them gently, and switch- 
ing on or off the effect of the side walls gave no qualitative 
difference. 

We found that both the angle of repose and the sur- 
face roughness exhibits linear dependence on the ratio / 
of the maximum magnetic force at contact and the grav- 
itational force. As it was also mentioned by Forsyth et 
al. P, 1 / overestimates the effective cohesion because 
of the anisotropy of the magnetic interaction. According 
to this, the angle of repose increases much more slowly 
with / than expected from stability criteria Q and ex- 
periments H on wet granular media. The experimental 
results of Forsyth et al. and our simulations are in good 
accordance, though the angle of repose at zero magnetiza- 
tion in our case was smaller. The side wall model taking 
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into consideration the effect of the front and back walls 
of a Hele-Shaw cell arrangement, could increase the sur- 
face angle but did also introduce a stronger dependence 
of the angle of repose on / than in the experiments of 
Forsyth et al. 0,0. Taking into consideration the static 
and rolling friction of particles could probably reproduce 
more closely the experimental results in this respect too. 

Tegzes et al. |5| found in experimental studies of wet 
granular media linear dependence of the angle of repose 
on the intcrparticle force ratio in the granular regime and 
an almost linear dependence in correlated regime with 
a slight curvature. This bending could not be clearly 
identified in our results and most probably would require 
more accurate investigations at both small and large in- 
terparticle force ratios. 

As reported by Tegzes et al. |6| in case of wet gran- 
ular media there is a difference in avalanche formation 
at small and at large interparticle force ratios. We could 
also identify a granular and a correlated regime in case 
of magnetized particles. The granular regime is charac- 
terized by quasi-continuous granular flows, while the cor- 
related regime is characterized by long particle clusters 
falling at the system's boundary. The transition between 
the two regimes is not sharp. In simulations we found 
that the transition is at f c ss 7, while calculations based 
on stability criteria indicate a transition at f c rs 6 for the 
investigated magnetic case. 

In the granular regime the avalanche sizes are propor- 
tional to the corresponding average avalanche durations. 
According to this, there is a well defined average granu- 
lar flow characterizing the avalanches. We found that this 
increases linearly with /. Analyzing the avalanche size 
distributions, we also found that there is a well-defined 
characteristic size in particle avalanches. Based on scal- 
ing properties of the avalanche size distributions, we ar- 
gue that the characteristic avalanche size increases lin- 
early with /. The characteristic average avalanche dura- 
tion in leading order seems to be independent of /. This 
dependence is contained in higher order corrections which 
could not be captured by our simulations. 

In the correlated regime the avalanche sizes are in lead- 



ing order proportional to the square of the corresponding 
average avalanche durations. This is explained by the free 
fall of long particle clusters. The avalanche durations are 
defined by the length of the falling particle clusters. The 
width of the falling clusters is small, and can take arbi- 
trary values up to a well defined maximum. According 
to this the avalanche sizes are defined in leading order by 
the cluster lengths, but there are large fluctuations which 
are also proportional to the cluster lengths. The cluster 
length increases linearly with /. The avalanche duration 
distributions show evidence of a characteristic avalanche 
duration, indicating also a characteristic avalanche size. 
We could not identify this avalanche size explicitly be- 
cause of the large fluctuations in the avalanche size dis- 
tribution. 

Our results regarding the avalanche size distributions 
in both granular and correlated regime are very close 
to the experimental results of Tegzes et al. t)j on wet 
granular materials. They could identify characteristic 
avalanche sizes in both regimes, and they state that the 
avalanche size distributions in correlated regime are broad 
and both small and large avalanches may occur. This 
seems to confirm our finding that in correlated regime 
there are large avalanche sizes fluctuations. 

The results on avalanche sizes and durations may 
slightly depend on the chosen time scale on which the 
avalanches are observed, however we argue that much 
more coarser or finer time scale will both lead to non- 
physical results, while small correction in the time scale 
will not lead to qualitative difference. 
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